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Abstract 


There are many costs associated with increased body size and longevity in animals, including the accumulation of genotoxic and 
cytotoxic damage that comes with having more cells and living longer. Yet, some species have overcome these barriers and have 
evolved remarkably large body sizes and long litespans, sometimes within a narrow window of evolutionary time. Here, we dem- 
onstrate through phylogenetic comparative analysis that multiple turtle lineages, including Galapagos giant tortoises, concurrently 
evolved large bodies, long lifespans, and reduced cancer risk. We also show through comparative genomic analysis that Galapagos 
giant tortoises have gene duplications related to longevity and tumor suppression. To examine the molecular basis underlying 
increased body size and lifespan in turtles, we treated cell lines from multiple species, including Galapagos giant tortoises, with 
drugs that induce different types of cytotoxic stress. Our results indicate that turtle cells, in general, are resistant to oxidative stress 
related to aging, whereas Galapagos giant tortoise cells, specifically, are sensitive to endoplasmic reticulum stress, which may give 
this species an ability to mitigate the effects of cellular stress associated with increased body size and longevity. 
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Significance 


Although all animals are sensitive to developing cancer, some are remarkably resistant to developing cancer during 
their lifetimes. Identifying how these species evolved to be cancer resistant can provide new insights into basic cancer 
biology and evolution. We found that long-lived turtles and tortoises evolved enhanced cancer resistance, at least in 
part by gaining extra copies of cancer-fighting genes. These observations may help the development of new cancer 
treatment and prevention strategies. 


Introduction 

Body size and longevity are fundamental lite history traits that blue whale (Balaenoptera musculus) (Magalhaes and Costa 
vary tremendously across vertebrates. Maximum body mass In 2009), whereas maximum lifespan ranges from 8 weeks in 
vertebrates ranges from 0.5g in the red-backed salamander the pygmy goby (Eviota sigillata) (Depczynski and Bellwood 


(Plethodon cinereus) (Moore et al. 2001) to 136,000 kg in the 2005) to over 400 years in the Greenland shark (Somniosus 
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microcephalus) (Nielsen et al. 2016). Life history comparisons 
also show a strong positive correlation between body size and 
lifesoan across animals, with few exceptions (Healy et al. 
2014). There are significant physiological constraints acting 
on organisms at the larger, longer-lived end of this spectrum, 
particularly the accumulation of mutations and cellular dam- 
age that comes with having more cells and greater cell turn- 
over (Peto 2015). The consequences of such long-term 
genotoxic and cytotoxic stress include genome instability, mi- 
tochondrial dystunction, telomere reduction, and increased 
cancer risk (LO6pez-Otin et al. 2013). 

A recurring theme in lifespan and aging regulation Is the 
critical role played by processes that promote cellular protec- 
tion and maintenance (Kenyon 2010), including the ability of 
cells to recycle materials, repair damage, and remove waste. 
Senescent cells, whose numbers greatly increase with age, 
exhibit declines in these processes, and are also associated 
with proinflammatory phenotypes that are linked to age- 
related diseases (Baar et al. 2017; Flatt and Partridge 2018). 
At the same time, apoptosis, which is the programmed de- 
struction of unfit or damaged cells, is reduced in older indi- 
viduals (Salminen et al. 2011). This decline in cell performance 
in combination with a decreased ability to remove poor- 
performing cells is central to the aging process (LOpez-Otin 
et al. 2013). Similarly, cancer can arise trom cumulative gen- 
otoxic and cytotoxic stress, and apoptosis also plays a primary 
role in cancer resistance by removing potentially cancerous 
cells (Verfaillie et al. 2013). Thus, if cancer-suppressing mech- 
anisms are similar across species, then larger, longer-lived 
organisms should be at greater risk of cancer than smaller, 
shorter-lived ones (Peto 2015). Although this correlation exists 
within species, for example, cancer incidence increases with 
increasing adult height for most cancer types in humans 
(Green et al. 2011) and overall body mass in dogs (Dobson 
2013), there is no such correlation between species—an ob- 
servation often referred to as “Peto’s paradox” (Peto 2015). 

The molecular and cellular mechanisms underlying the evo- 
lution of large bodies and long lifesoans have been explored in 
mammals such as elephants (Sulak et al. 2016; Vazquez et al. 
2018), whales (Seim et al. 2014), bats (Foley et al. 2018; 
Gorbunova et al. 2020), and naked mole rats (Salmon et al. 
2008; Gorbunova et al. 2014), but are less well studied in 
other vertebrates. Reptiles are an excellent system in which 
to study the evolution of body size and longevity because 
diverse lineages have repeatedly evolved large body sizes 
and long lifespans (Chiari et al. 2018). Turtles, in particular, 
have lower rates of neoplasia than snakes and lizards (Garner 
et al. 2004; Sykes and Trupkiewicz 2006), are especially long- 
lived, and are “slower aging” than other reptiles (Hoekstra 
et al. 2020). Most notably, Galapagos giant tortoises 
(Chelonoidis niger species complex; hereafter referred to as 
C. niger) and Aldabra giant tortoises (A/dabrachelys gigantea) 
can live over 150 years (3-5 times longer than their closest 
relatives) and weigh over 200 kg (50-100 times heavier than 
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their closest relatives) (Caccone et al. 1999; Palkovacs et al. 
2002; Poulakakis et al. 2012; Chiari 2020). Galapagos giant 
tortoises also appear to have evolved a suite of cellular traits 
that may contribute to their longevity, such as a slower rate of 
telomere shortening and extended cellular litesoans com- 
pared with mammals (Goldstein 1974). 

Here, we explore the evolution of body size and lifespan in 
turtles by integrating several approaches (fig. 1): 1) phyloge- 
netic comparative analysis of body size, lifespan, and intrinsic 
cancer risk in turtles; 2) gene duplication analysis of aging and 
cancer-related genes across available turtle genomes; and 3) 
cell-based assays of apoptosis and necrosis in multiple turtle 
species varying in body size and lifespan. We show that spe- 
cies with remarkably long lifespans, such as Galapagos giant 
tortoises, also evolved reduced cancer risk. We also confirm 
that the Galapagos giant and desert tortoise genomes encode 
numerous duplicated genes with tumor suppressor and anti- 
aging functions (Quesada et al. 2019). Our comparative ge- 
nomic analysis further suggests that cells trom large, long- 
lived species may respond differently to cytotoxic stress, in- 
cluding endoplasmic reticulum (ER) and oxidative stress. The 
combined genomic and cellular results suggest that at least 
some turtle lineages evolved large bodies and long lifespans, 
in part, by increasing the copy number of tumor suppressors 
and other antiaging genes and undergoing changes in cellular 
phenotypes associated with cellular stress. 


Results 


Repeated Evolution of Large Body Size in Turtles 


We found substantial independent accelerations in the rate of 
body size evolution in several turtle lineages, including a 29x 
rate increase (386% increase in carapace length) in the stem- 
lineage of sea turtles (Chelonioidea) and a further 103~x rate 
increase (757% increase in carapace length) in leatherback 
sea turtles (Dermochelys coriacea), a 37x rate increase (35% 
increase in carapace length) in the stem-lineage of soft-shell 
turtles (Trionychidae), a 200 x rate increase (364% increase in 
carapace length) in the stem-lineage of narrow-headed soft- 
shell turtles (Chitra chitra and Chitra indica), and a 463 rate 
increase (364% increase in carapace length) in Cantor's giant 
softshell turtle (Pelochelys cantoril) (fig. 2A). 

Among the more notable groups with increased rates of 
body size evolution were the “giant” tortoises, including a 
81x rate increase (10% increase in carapace length) in the 
stem-lineage of recently extinct Mascarene giant tortoises 
(Cylindraspis spp.), a 137x rate increase (55% increase in 
carapace length) in the stem-lineage of Aldabra giant tortoises 
(Aldabrachelys spp.) and a 87x rate increase (383% increase 
in carapace length) in A. gigantea, and a series of rate accel- 
erations in the ancestral lineages of Galapagos giant tortoises, 
including a 20x rate increase (19% increase in carapace 
length) in the stem-lineage of Geochelone and Chelonoidis, 
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Fic. 1.—Overview of the study design. Species with genomes utilized for gene duplication analysis as well as turtle species with cells used to measure 
apoptotic responses to genotoxic and cytotoxic drugs are indicated. Phylogenetic tree was built with TimeTree (Kumar et al. 2017). Turtle size data come 
from Ernst and Barbour (1992) and Colston et al. (2020), longevity data are from AnAge. 


a 20x rate increase (27% increase in carapace length) in the 
stem-lineage of Chelonoidis, and a 49x rate increase (81% 
increase in carapace length) in the stem-lineage of C. niger. 
These data indicate that gigantism evolved independently in 
multiple lineages of turtles, and step-wise in the evolution of 
giant tortoises with several rate accelerations in lineages an- 
cestral to C. niger. 


Reduction of Intrinsic Cancer Risk in Turtles 


In order to account for a relatively constant prevalence of can- 
cer across species (Dorn et al. 1968; Abegglen et al. 2015; 
Boddy et al. 2020), intrinsic cancer risk must coevolve with 
changes in body size and lifespan across species. For example, 
a 100-year retrospective study of neoplasia in zoo reptiles iden- 
tified only six neoplasms in 490 turtle necropsies, which ranged 
in size from the West African mud turtle (Pelusios castaneus, 
carapace length ~25-28cm) to the spiny softshell turtle 
(Apalone spinifer spinifer, carapace length ~54cm) (Sykes 
and Trupkiewicz 2006). As expected, relative intrinsic cancer 
risk (RICR) in turtles also varies with changes in body size and 
lifespan (Tig. 2A). We estimated a 73-logz decreased RICR in 
the stem-lineage of sea turtles (Chelonioidea), a 129-log> de- 
creased RICR in leatherback sea turtles (Dermochelys coriacea), 
a 14-log> decrease in RICR in the stem-lineage of soft-shell 
turtles (Trionychidae), a 34-logz decreased RICR in the stem- 
lineage of narrow-headed softshell turtles (Chitra chitra and 
Chitra indica), and a 140-logz decreased RICR in the stem- 
lineage of Cantor’s giant softshell turtle (Pelochelys cantoril). 
Among the “giant” tortoises, we estimated a 97-log> de- 
creased RICR in the stem-lineage of Mascarene giant tortoises 
(Cylindraspis spp.), a 154-logy decreased RICR in the stem- 
lineage of Aldabra giant tortoises (A/dabrachelys spp.) and a 
50-logz decreased RICR in A. gigantea. In the lineages ances- 
tral to Galapagos giant tortoises, we estimated a 67-log> 
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decreased RICR in the stem-lineage of C. niger a 27-logz de- 
crease RICR in the stem-lineage of the Chelonoidis, and a 19- 
log decreased RICR in the stem-lineage of Geochelone and 
Chelonoidis. Thus, turtles coevolved large bodies and reduced 
intrinsic cancer risk, including step-wise reductions in the lin- 
eages ancestral to C. niger. 


Identification of Tumor Suppressor and Antiaging Gene 
Duplications in Turtle Genomes 


Previous studies have shown that large-bodied cancer-resistant 
species such as elephants (Sulak et al. 2016; Vazquez and Lynch 
2021) and whales (Keane et al. 2015) evolved an increased 
number of tumor suppressors, suggesting that the same may 
be possible in giant, long-lived turtles. A previous study of 
Galapagos giant tortoises, for example, identified several gene 
duplications in pathways that might be related to body size 
evolution and reduced cancer risk (Quesada et al. 2019). 
Therefore, we reanalyzed the Galapagos giant tortoise genome 
and other turtle genomes to Identify gene duplications and 
used maximum likelinood-based ancestral state reconstruction 
(ASR) to determine lineages in which genes were duplicated. 
We identified ~86 duplications in the stem-lineage of 
Testudines, 245 in the stem-lineage of Pleurodira, 33 in the 
stem-lineage of tortoises (AncTortoise), 259 in Chelonoidis 
abingdonii, 201 in Gopherus agassizii, 273 in Terrapene caro- 
lina, 315 in C picta, and 270 in Pelodiscus sinensis (fig. 3A). 
Consistent with previous studies which observed duplica- 
tion of tumor suppressor and other antiaging genes in large, 
long-lived species, we found that 12% of the pathways 
enriched among Galapagos giant tortoises were related to 
cancer and aging biology, whereas only 0O-6% of the path- 
ways that were enriched among gene duplications in other 
lineages were related to cancer and aging biology (fig. 38). 
Next, we identified GO cellular component terms that were 
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Fic. 2.—Convergent evolution of large-bodied, cancer-resistant turtles. (A) Turtle phylogeny with branch lengths and colors scaled by logz change in 
carapace length as a proxy for body size (left) and estimated intrinsic cancer risk (right). Clades and lineages leading to exceptionally large turtles and tortoises 
are labeled, as are species used in cytotoxic stress assays. (B) Lifespan (upper), body size (middle), and estimated intrinsic cancer risk (lower) of species used in 


stress assays. 


enriched among gene duplications in each lineage. Although 
gene duplications in some lineages were enriched in GO 
terms related to cancer biology and aging, significantly 
more ontology terms in Galapagos giant tortoises were re- 
lated to cancer and aging biology (fig. 3C). Enriched pathway 
and ontology terms in Galapagos (fig. 3D) tortoises included 
“apoptosis,” “programmed cell death,” “cell death signal- 
ling via NRAGE, NRIF and NADE,” “dual incision in GG-NER” 
and “formation of incision complex in GG-NER,” and 
“reduction of cytosolic Ca++ levels.” We also observed 


W i 


that “regulation of insulin-like growth factor (IGF) transport 
and uptake by Insulin-like Growth Factor Binding Proteins 
(IGFBPs)” was an enriched pathway term among 
Galapagos giant tortoise gene duplications, which may be 
related to the regulation of body size. Among the GO cellular 
component terms exclusively enriched among Galapagos gi- 
ant tortoise gene duplications were “ER membrane protein 
complex,” “endoplasmic reticulum membrane” and 
“nuclear outer membrane-endoplasmic reticulum membrane 
network,” and “anaphase-promoting complex.” 
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Fic. 3.—Gene duplicates in Galapagos and desert tortoises are enriched in tumor suppressor and antiaging functions. (A) Turtle phylogeny indicating the 
number of genes duplicated in each lineage, inferred by maximum likelihood. Inset, phylogeny with branch lengths proportional to gene duplication rate. 
The asterisk (*) denotes a node with gene duplications reconstructed with lower support than other nodes and nonsignificant (BPP=0.541). (B) Pie charts 
indicating the proportion of enriched Reactome pathways in each lineage related to cancer biology and aging (blue slices). Gene duplicates in Galapagos 
giant and desert tortoises are significantly more enriched in these terms than other lineages (two-sided permutation t-test is 0.00). (C) Pie charts indicating the 
proportion of enriched GO cellular component terms related to cancer biology, DNA damage repair, programmed cell death, and the endoplasmic reticulum 
(red slices). Gene duplicates in Galapagos giant and desert tortoises are significantly more enriched in these terms than other lineages (two-sided permutation 
t-test is 0.00). (D) Wordcloud of the Reactome (blue) pathways and GO cellular component terms (red) enriched exclusively in Galapagos giant tortoises. Only 
pathway and GO terms enriched with P<0.05 are shown are scaled according to log>-fold enrichment (see inset scale). (F) Wordcloud of the Reactome (blue) 
pathways and GO cellular component terms (red) enriched in desert tortoises. Only pathway and GO terms enriched with P<0.05 are shown are scaled 


according to log2-fold enrichment (see inset scale). 


Desert tortoise-specitic gene duplications were enriched in 
pathways related to cancer biology and aging, particularly 
compared with other turtles. For example, 9.4% of the path- 
ways enriched among desert tortoise duplicates were related 
to cancer and aging biology, significantly more than other tur- 
tle lineages but less than Galapagos giant tortoises (fig. 38). 
Similarly, significantly more GO terms in desert tortoises were 
related to cancer and aging biology (fig. 3C). Enriched pathway 
and ontology terms in desert tortoises (fig. 3£) were related to 
DNA damage and repair including “Formation of TC-NER Pre- 
Incision Complex,” “ Gap-tilling DNA repair synthesis and liga- 
tion in TC-NER,” “Dual incision in TC-NER,” “Apoptotic cleav- 
age of cell adhesion proteins,” and “Transcription-Coupled 
Nucleotide Excision Repair (TC-NER).” Enriched in GO terms 
included “integral component of mitochondrial outer mem- 
brane,” = “nucleotide-excision repair complex,” and 
“autophagosome.” These data suggest that desert tortoises 
have evolved gene duplications that may also contribute to 
cancer resistance and the evolution of longevity. 
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Turtle Cells Have Unique Responses to Genotoxic and 
Cytotoxic Stress 


Our observation that gene duplications in the Galapagos 
and desert tortoise genomes are enriched in pathways and 
GO terms related to the biology of aging, apoptosis, cell 
cycle regulation, DNA damage repair, and mitochondrial ox- 
idative DNA damage protection (fig. 3D), suggests that cells 
trom these species may have different cellular responses to 
genotoxic and cytotoxic stress than cells from other turtles. 
To test this hypothesis, we treated primary fibroblasts trom 
C. niger, Geochelone platynota, Gop. agassizii, Homopus 
aerolatus, and T. carolina (fig. 2B) with drugs to induce dif- 
ferent types of stress including: 1) tunicamycin, which indu- 
ces ER stress and the unfolded protein response (UPR) 
through an accumulation of unfolded and misfolded pro- 
teins (Banerjee et al. 2011; Guha et al. 2017); 2) etoposide, 
which forms a ternary complex with DNA and topoisomer- 
ase Il and prevents religation of replicating DNA strands 
leading to single- and double-stranded DNA _ breaks 
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Fic. 4.—Cells from Galapagos giant and desert tortoises have unique stress responses. (A) Left, dose response curves for tunicamycin, which induces 
endoplasmic reticulum stress and the unfolded protein response (24 h posttreatment). Right, boxplots showing differences between species 24 h after 
treatment with 25 wM tunicamycin; statistical tests are relative to C. niger, which had the strongest apoptotic response. The unpaired mean difference 
between C. niger and G. platynota is —0.521 (95% Cl —0.631, —0.408). The P value of the two-sided permutation t-test is 0.0. The unpaired mean 
difference between C. niger and Gop. agassizii is —0.288 (95% Cl —0.399, —0.193). The P value of the two-sided permutation t-test is 0.0. The unpaired 
mean difference between C. niger and H. areolatus is —0.278 (95% Cl —0.431, —0.124). The P value of the two-sided permutation t-test is 0.0. The 
unpaired mean difference between C. niger and T. carolina is —0.355 (95% Cl —0.465, —0.246). The Pvalue of the two-sided permutation t-test is 0.0. n=3. 
(B) Left, dose response curves for etoposide, which induces DNA strand breaks (24 h posttreatment). Right, boxplots showing differences between species 
24 h after treatment with 500 uM etoposide; statistical tests are relative to Gop. agassizii, which had the strongest apoptotic response. The unpaired mean 
difference between Gop. agassizii and C. niger is —0.141 (95% Cl —0.195, —0.106). The P value of the two-sided permutation t-test is 0.0. The unpaired 
mean difference between Gop. agassizii and G. platynota is —0.16 (95% Cl —0.227, —0.112). The P value of the two-sided permutation t-test is 0.0. The 
unpaired mean difference between Gop. agassizii and H. areolatus is —0.16 (95% Cl —0.227, —0.112). The P value of the two-sided permutation t-test is 
0.0. The unpaired mean difference between Gop. agassizii and T. carolina is —0.147 (95% Cl —0.198, —0.115). The P value of the two-sided permutation t- 
test is 0.0. n=3. (C) Left, dose response curves for paraquat, which induces oxidative stress (24 h posttreatment). Right, boxplots showing differences 
between species 24 h after treatment with 2,500 uM paraquat; statistical tests are relative to C. niger, which had the weakest apoptotic response. The 
unpaired mean difference between C. niger and G. platynota is 0.24 (95% Cl 0.176, 0.317). The P value of the two-sided permutation t-test is 0.0. The 
unpaired mean difference between C. niger and Gop. agassizii is 0.272 (95% Cl 0.188, 0.356). The P value of the two-sided permutation t-test is 0.0. The 
unpaired mean difference between C. niger and H. areolatus is 0.0927 (95% Cl 0.0337, 0.169). The P value of the two-sided permutation t-test is 0.0. The 
unpaired mean difference between C. niger and T. carolina is 0.32 (95% Cl 0.262, 0.396). The P value of the two-sided permutation t-test is 0.0. n=3. 


(Wozniak and Ross 1983); and 3) paraquat, which causes We found that tunicamycin induced a dose-dependent in- 
oxidative stress through the production of reactive oxygen crease in apoptosis in cells trom most species (fig. 4A). 
species and S-phase cell cycle arrest (Salmon et al. 2008). However, at 24h and 25 uM, C niger cells had an apoptotic 
We quantified the kinetics of cell death using the RealTime- response that was at least double that of other species, 
Glo™ Annexin V Apoptosis and Necrosis assay (RTG) every whereas G. platynota cells were insensitive to tunicamycin 
30 min for 48h. (fig. 4A). In contrast etoposide did not induce apoptosis or 
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necrosis in cells trom most species but did induce a strong 
apoptotic response in Gop. agassizii cells (fig. 48). Turtle cells 
were also variably sensitive to paraquat, but cells trom all 
species induced apoptosis in response to paraquat treatment 
(fig. 4C). Thus, C. niger cells are more sensitive to ER and UPR 
stress than other species, Gop. agassizii cells are more sensi- 
tive to DNA damage induced by etoposide than other species, 
and cells from all species were sensitive to oxidative stress 
induced by paraquat. 


Discussion 


Reptiles in general, and turtles specifically, are an excellent 
system for studying the mechanisms underlying variation in 
body size, lifesoan, and cancer resistance (Chiari et al. 2018; 
Hoekstra et al. 2020). Adult body size among turtle species 
differs by three orders of magnitude, ranging from the 
speckled dwarf tortoise (Chersobius signatus; 100g) to 
Aldabra giant tortoises (A. gigantea; >300 kg), whereas lon- 
gevity across species varies from tens of years to > 150 years in 
giant tortoises. Yet, even the smallest turtles live relatively long 
(up to 30years) compared with most other vertebrates. 
Turtles also have lower estimated cancer rates (~ 1.2%) com- 
pared with mammals (~12.5%), suggesting that they have 
evolved the means to delay aging and reduce cancer suscep- 
tibility (Garner et al. 2004; Sykes and Trupkiewicz 2006; 
Abegglen et al. 2015; Boddy et al. 2020). Although increased 
lifespan and cancer prevalence in turtles could, in part, be due 
to reduced damage resulting from their lower metabolic rates, 
previous genomic and cellular data suggest that there may 
also be molecular differences that enable extremes in body 
size, longevity, and cancer resistance (Goldstein 1974: 
Quesada et al. 2019). 

Here, we found that body size rapidly increased in multiple 
turtle lineages independently, including a dramatic increase In 
both body size and the rate of body size evolution in 
Galapagos giant tortoises. When body size, longevity, and 
intrinsic cancer rates are analyzed together within an explicit 
phylogenetic context, large-bodied lineages evolved reduced 
intrinsic cancer risk. Among the most dramatic decreases in 
cancer risk among turtles was in the stem-lineage of 
Galapagos giant tortoises, suggesting this lineage evolved ge- 
netic and cellular mechanisms that reduce cancer risk. 
Consistent with the evolution of reduced intrinsic cancer risk 
in this lineage, we found that gene duplications in the 
Galapagos lineage are enriched in ER stress associated path- 
ways Is similar to a previous analysis by Quesada et al. (2019), 
who found that genes under positive selection in Galapagos 
giant tortoises are also enriched for ER function and stress 
pathways. ER stress dysregulates protein homeostasis, which 
is one of the hallmarks of aging and reduced animal lifespan 
(LOpez-Otin et al. 2013; Morimoto and Cuervo 2014; Kaushik 
and Cuervo 2015). The impairment of protein production, 
folding, and degradation associated with ER stress can impact 
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many cellular processes and lead to the buildup of toxic pro- 
tein aggregates that are linked to many age-related diseases 
and cancers (Vilchez et al. 2014). Our observations that du- 
plicate genes are enriched in anticancer and antiaging path- 
ways Is similar to several previous studies (Caulin et al. 2015; 
Tollis et al. 2020; Tejada-Martinez et al. 2021; Vazquez and 
Lynch 2021), suggesting that increased copy number of these 
genes may be a general trend in large, long-lived species. 

Our data suggest a mechanistic connection between genes 
related to the ER, enhanced responses to ER stress, and the 
UPR in Galapagos giant tortoise cells. Indeed, we found that 
tunicamycin, which induces ER stress by activating the UPR 
and ultimately inducing apoptosis (Nami et al. 2016), causes 
an immediate and pronounced apoptotic response in 
Galapagos giant tortoise cells, whereas cells from other spe- 
cies were much slower to respond. These results indicate that 
Galapagos giant tortoise cells have an extremely sensitive ER 
stress response. In contrast, Galapagos giant tortoise cells did 
not exhibit a heightened apoptotic response to paraquat and 
etoposide, which induce oxidative stress, DNA damage, and 
apoptosis independently of the ER stress/UPR-signaling path- 
way (Mizumoto et al. 1994; Suntres 2002). These data sug- 
gest that changes In the ER stress and UPR-signaling pathways 
may contribute to the evolution of long lifespans, large bod- 
ies, and augmented cancer resistance in Galapagos giant 
tortoises. 

Unexpectedly, we also found that cells from most species 
were unresponsive to etoposide treatment, even at the high- 
est dose and longest exposure time. The sole exception was 
Gop. agassizii cells, which were similar to other species in their 
response to tunicamycin and paraquat, but markedly more 
sensitive to etoposide treatment. These data suggest that tur- 
tle cells can either export intracellular etoposide before It can 
induce DNA damage, are generally insensitive to DNA dam- 
age induced by etoposide, or can rapidly repair etoposide- 
induced DNA damage. Regardless of these, and potentially 
other mechanisms that alter sensitivity to etoposide, Gop. 
agassizii cells resoond differently than cells trom the other 
species tested. The functional and organismal consequences 
of this altered sensitivity are unclear, but is similar to elephant 
cells which have evolved to induce apoptosis at relatively low 
levels of DNA damage(Sulak et al. 2016; Vazquez et al. 2018; 
Vazquez and Lynch 2021). This reduced threshold for induc- 
iNg apoptosis in response to cellular stress may clear cells that 
have been exposed to the kinds of stresses that eventually 
lead to cancer before transformation into cancer cells. 

Although medical literature has often portrayed high levels 
of apoptosis as a maladaptive response to cellular stress con- 
tributing to aging (Szegezdi et al. 2006; Chadwick and Lajoie 
2019), previous work in elephants showed that heightened 
apoptotic responses to genetic and cellular damage can ac- 
tually be an adaptive and beneficial response linked to in- 
creased body size and longevity (Abegglen et al. 2015; 
Sulak et al. 2016; Vazquez et al. 2018). This is because rapid 
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and effective clearance of damaged or injured cells can help 
maintain tissue integrity (Baar et al. 2017; de Keizer 2017), 
especially in organisms with many cells and large amounts of 
cell turnover. Thus, our results are compatible with ER stress as 
a potential factor in the evolution of large, long-lived turtles. 


Caveats and Limitations 


There are several limitations to our findings that could be the 
subject of Tuture research. Although all cellular assays were 
performed on a common cell type, fibroblasts, the tissue or- 
igin of these cells differed among species. We note, however, 
that for Galapagos giant tortoises, tunicamycin response was 
not confounded with the site of fibroblast origin. We were 
also unable to obtain the age of source animals, which could 
affect results since apoptosis can decline during senescence 
(Salminen et al. 2011). We also have not classified duplicate 
genes into functional categories, such as those that have pre- 
served ancestral functions, evolved new ones, or evolved new 
expression domains because these classifications are ex- 
tremely difficult to predict based on sequence characteristics 
alone. Finally, the taxon sampling can be expanded to include 
cells from closely related species with large differences in body 
size or lifespan. This would enable better resolution in isolat- 
ing the evolutionary origins of enhanced responses to genetic 
and cellular stress, for example, by including the closest living 
relative of Galapagos giant tortoises (Chaco tortoises; 
Chelonoidis chilensis). We relied on primary cell lines that 
were currently available in frozen zoos and commercial bio- 
banks; but future cellular work could attempt to sample fresh 
cells trom the same tissue type and life stage across species, 
which is logistically challenging. 


Conclusions 


Although numerous studies have found comparative genomic 
signatures associated with the evolution of body size, longev- 
ity, and cancer resistance (Keane et al. 2015; Herrera-Alvarez 
et al. 2021; Babarinde and Saitou 2020), there have been few 
attempts to validate these findings experimentally at the cel- 
lular level Vimenez et al. 2018). Furthermore, most previous 
cellular studies on these subjects focus almost exclusively on 
placental mammals. The work presented here utilizes turtle 
cell lines, which is much more feasible than studying body size 
and aging phenotypes in these long-lived animals. Our most 
salient finding is that Galapagos giant tortoises are much 
more sensitive at inducing apoptosis in response to ER stress 
compared with other turtle species, and also have genomic 
and phylogenetic signatures of rapid evolutionary increases in 
body size, lifesoan, and cancer resistance. We also found 
more generally that all turtle cell lines were resistant to oxida- 
tive stress induced by paraquat. This supports previous oxida- 
tive stress studies in turtles, and indicates that turtles, in 
general, may be a promising model system in which to study 


GBE 


resistance to stress trom long lifespans (Lutz et al. 2003). 
However, although the resistance of turtle cells to oxidative 
stress may contribute to the generally long lifespans of turtles, 
the sensitivity of Galapagos giant tortoise cells to ER stress 
may increase their resistance to oncogenic transformation 
thus promoting healthy aging and larger body sizes. 


Materials and Methods 


Intrinsic Cancer Risk Estimation 


The dramatic increase in body size and lifespan in some turtle 
lineages, and the relatively constant rate of cancer across spe- 
cies of diverse body sizes and lifespan (Leroi et al. 2003), 
would predict an increase in cancer risk concurrent with an 
increase in body size or lifespan. In order to identity lineages 
with exceptional changes in body size (with total carapace 
length used here as a proxy for body size), longevity, or in- 
trinsic cancer risk, we jointly estimated these parameters 
across turtles and reconstructed ancestral states within a phy- 
logenetic framework. Using phylogenetic and body size data 
trom Colston et al. (2020), which includes 357 extant turtle 
species five species of extinct Cylindraspis, and longevity data 
trom AnAge (Magalhaes and Costa 2009), and following 
Peto’s (2015) model of cancer risk (Vazquez and Lynch 
2021), we estimated the intrinsic cancer risk (K) as the product 
of risk associated with body size and lifespan (lifespan°x body 
size). In order to determine (K) across species and at ancestral 
nodes, we first estimated body size at each node. We used a 
generalization of the Brownian motion model that relaxes 
assumptions of neutrality and gradualism by considering 
increments to evolving characters to be drawn from a 
heavy-tailed stable distribution (the stable model) imple- 
mented in StableTraits (Elliot and Mooers 2014). The stable 
model allows for large jumps in traits and has previously been 
shown to out-perform other models of body size evolution, 
including standard Brownian motion models, Ornstein— 
Uhlenbeck models, early burst maximum likelihood models, 
and heterogeneous multirate models (Prang 2019). We used 
phylogenetic generalized least-square regression (Grafen and 
Hamilton 1989; Martins and Hansen 1997; Pagel 1997) using 
a Brownian covariance matrix as implemented in the R pack- 
age ape (Paradis and Schliep 2019) to inter ancestral lifespans 
across turtles using our estimates for body size (Colston et al. 
2020) and reported maximum lifespans for each species 
(Magalhaes and Costa 2009). Fold-change in cancer suscep- 
tibility was estimated between all ancestral (K1) and descen- 
dant (K2) nodes. The fold-change in cancer risk between a 
node and its ancestor was then defined as K2/K1 (supplemen- 
tary fig. 2—source data 1, Supplementary Material online). 
We note that two variables that are correlated with each 
other, such as body mass and lifespan, will be correlated with 
the product of those variables, such as intrinsic cancer risk. 
Thus, we may have reached similar conclusions about 
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changes in interred cancer risk across turtle evolution using 
either body size or lifesoan data alone. But, we would have 
not inferred the magnitudes of those changes correctly for 
most species and may have inferred the wrong direction for 
some species that decrease dramatically in body size but in- 
crease dramatically in lifespan (although there are few such 
species). 


Identification of Duplicated Genes and Reconstruction of 
Ancestral Copy Numbers 


Following Caulin (2015), we identified duplicated genes in the 
genomes of the Pinta Island Galapagos giant tortoise (C. 
abingdonir; ASM359739v1), Goode’s thornscrub tortoise 
(Gopherus evgoodel; rGopEvg1_v1.p) and Agassiz’s desert 
tortoise (Gop. agassizil; ASM289641v1), Chinese softshell tur- 
tle (Pelodiscus sinensis; PelSin_1.0), painted turtle (Chrysemys 
picta belli, Chrysemys_picta_bellii-3.0.3), three-toed box tur- 
tle (T. carolina triunguis; T_m_triunguis-2.0), chicken (Gallus 
gallus; GRCg6a), and eastern brown snake (Pseudonaja texti- 
lis; EBS10Xv2-PRI) using the Ensembl (Genes 103) BioMart 
web-based tool to extract same-species paralogies (within_s- 
pecies_paralog) from each genome; same-species paralogies 
are identified by Ensembl Compara (Howe et al. 2021) using 
gene tree species tree reconciliation. Specifics of the Ensembl 
Compara (Howe et al. 2021) method can be found in the 
Ensembl Compara documentation: https://useast.ensembl. 
org/info/genome/compara/index.html. Genes that were clas- 
sified as pseudogenes on Ensembl were not included. 
Duplicate genes identified in Goode’s thornscrub tortoise 
(Gop. evgoode/; rGopEvg1_v1.p) were manually verified in 
Agassiz’s desert tortoise (Gop. agassizii; ASM289641v1) using 
reciprocal best BLAT using default parameters: https://useast. 
ensembl.org/Help/View?id=45 1#toc6. Interested readers are 
referred to these documents for specifics about the Compara 
and BLAST/BLAT methods as a complete description is beyond 
the scope of this manuscript. 

We used maximum likelinood-based ASR to determine 
when in the evolution of turtles each gene was duplicated. 
We encoded the copy number of each putatively functional 
gene for each species as a discrete trait, with state O for one 
gene copy and state 1 tor two or more copies. We used IQ- 
TREE to select the best-fitting model of character evolution 
(Kalyaanamoorthy et al. 2017; Minh et al. 2020), which was 
inferred to be a general time reversible model for morpholog- 
ical data (GTR2) with character state frequency optimized (FO) 
by maximum-likelihood from the data. Next, we inferred gene 
duplication events with the empirical Bayesian ASR method 
implemented in IQ-TREE (Kalyaanamoorthy et al. 2017; Minh 
et al. 2020), the best-fitting model of character evolution 
(GTR2+FO), and the unrooted species tree for turtles 
(fig. 3A). We considered ASRs to be reliable if they had 
Bayesian posterior probability (BPP) greater than or equal to 
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observed state frequency from the alignment; less reliable 
reconstructions were excluded from further analyses. 


Gene Duplication Pathway Enrichment Analysis 


To determine if gene duplications were enriched in particular 
biological pathways, we used Enrichr (Chen et al. 2013, p. 5; 
Kuleshov et al. 2016) to perform overrepresentation analysis 
(ORA) of the Reactome database. Furthermore, we used the 
Panther GO enrichment analysis tool (Mi et al. 2021) to per- 
torm ORA on GO cellular component terms (Ashburner et al. 
2000; Gene Ontology Consortium 2021). Gene duplicates in 
each lineage were used as the foreground gene set, and the 
initial query set was used as the background gene set. Enrichr 
uses a hypergeometric test for statistical significance of path- 
way overrepresentation, whereas the Panther GO enrichment 
analysis tool uses a binomial test for statistical significance of 
GO cellular component term overrepresentation. 


Cell Culture 


Experimental cellular phenotypes were generated to compare 
results from the body size, longevity, and intrinsic cancer risk 
analysis and the gene duplication and enrichment analysis. 
Apoptosis was chosen as the major cellular endpoint of inter- 
est because of its central role in aging and cancer through the 
removal of damaged or cancerous cells (Salminen et al. 2011; 
Verfaillie et al. 2013). 

We cultured cell lines trom four tortoise taxa (C. niger, G. 
platynota, Gop. agassizii, H. aerolatus) obtained from the San 
Diego Frozen Zoo and one turtle species (7. carolina) obtained 
trom the American Type Culture Collection (table 1). All cells 
were primary fibroblasts derived from either the heart 
(T. carolina), trachea (Gop. agassizii, G. platynota, C. niger), 
or eye (H. aero/atus). Turtle cells have previously been shown 
to grow within a range of 23-30 °C (Clark and Karzon 1967; 
Clark et al. 1970; Goldstein 1974). Therefore, cells were in- 
cubated at 25 °C, with 5% COs. Cells were cultured in min- 
imum essential medium (MEM; Gibco) with 10% fetal bovine 
serum (Gibco) and 1% _ penicillin-streptomycin antibiotic 
(Gibco) in standard T75 flasks (Thermo Fisher Scientific). 
Media was changed every 3 days. Cells were passaged before 
reaching 90% confluency, approximately every 7—9 days. For 
passaging, cell plates were rinsed with one volume of 37 °C 
DPBS (Gibco) and cells detached with 0.25% trypsin-EDTA 
(Gibco). We note that cells detached quickly with the assis- 
tance of gentle tapping of plates and without incubation. The 
cell suspension was transferred to a 15-ml conical tube 
(Thermo Fisher Scientific) with an equal volume of complete 
media to stop trypsinization. Cells were then centrifuged at 
500g for 5 min, and then the pellet was resuspended in 1 ml 
of complete media. Cell viability was determined using a 
TC10 Automated Cell Counter (Bio-Rad Laboratories) and 
was >75% for all cell lines throughout the experiment. 
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Table 1 
Turtle Experimental Cell Lines 
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Scientific Name Common Name Family Biopsy Site Cell Viability (%) 
Homopus areolatus Parrot-beaked tortoise Testudinidae Eye 83+ 15 
Gopherus agassizii Desert tortoise Testudinidae Trachea 93+6 
Geochelone platynota Burmese star tortoise Testudinidae Trachea 96+ 1 
Chelonoidis nigra Galapagos tortoise Testudinidae Trachea 89+7 
Terrapene Carolina Common box turtle Emydidae Heart 91+7 


Note.—Biopsy site is the location from which primary fibroblast cells were derived. Average in vitro cell viability percent is over 12-14 passages depending on cell line. Cell 


viability was calculated as the ratio of live cells to total cells at each passage with SD. 


We selected various cytotoxic drugs to induce different 
types of cellular stress, including etoposide (Cayman 
Chemical Company), which induces single- and double- 
stranded DNA breaks (Wozniak and Ross 1983), paraquat 
(Sigma-Aldrich), which induces oxidative stress through pro- 
duction of reactive oxygen species and s-phase cell cycle arrest 
(Salmon et al. 2008), and tunicamycin (Cayman Chemical 
Company), which induces ER stress and the UPR by causing 
an accumulation of unfolded and misfolded proteins 
(Banerjee et al. 2011; Guha et al. 2017). 


Kinetic Measurements of Cell Death 


The RTG visualizes the kinetics of apoptosis over a given pe- 
riod of time and differentiates secondary necrosis occurring 
during late apoptosis from necrosis caused by other cytotoxic 
events (Landreman et al. 2019). RTGs were performed for 
each species by seeding 5,000 cells per well into an opaque 
bottomed 96-well plate with three replicates per treatment 
per species and one empty column with no cells (background 
control). Cells were left to adhere for 24h, after which seed- 
ing media was aspirated off and serial dilutions of 
Tunicamycin (0, 10, 50, 100 uM), Etoposide (0, 100, 500, 
1,000 uM), or Paraquat (0, 100, 1,000, 5,000 wM) were ap- 
plied. All drug treatments were made using Fluorobrite 
DMEM media. The “0 uM” control treatments consisted of 
the vehicle used (DMSO or PBS) at the concentration match- 
ing the highest drug concentration, whereas the background 
control (NoCell) consisted of the drug treatment and assay 
reagents with no cells. As per the RTG product protocol, the 
500-fold dilution of reagents in Fluorobrite DMEM media was 
added to wells immediately after drug treatments were ap- 
plied (effectively reducing the initial drug dilutions in halt). 
Readings were then taken every 30min for 48h using a 
GloMax Luminometer (Promega Corporation). 


Statistical Analysis 


All statistical analyses reported in this paper are estimation 
Statistics, including effect sizes, 95% confidence intervals 
(Cls) of the effect size, and P values. Effect sizes and 95% 
Cls are reported as: effect size (CI width lower bound; upper 
bound) with 5,000 bootstrap samples; the Cl is bias-corrected 


and accelerated. P values reported are the likelihood of ob- 
serving the effect sizes if the null hypothesis of zero difference 
is true. For each permutation P value, 5,000 reshuffles of the 
control and test labels were performed (Ho et al. 2019). All 
Statistical analyses were performed in RStudio (RStudio Team 
2020). 


Supplementary Material 


Supplementary data are available at Genome Biology and 
Evolution online. 
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